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ABSTRACT 

We present equations for the gravitational lensing flexion expected for an elliptical 
lens mass distribution. These can be reduced to one-dimensional finite integrals, thus 
saving significant computing time over a full two-dimensional calculation. We estimate 
constraints on galaxy halo ellipticities for a range of potential future surveys, finding 
that the constraints from the two different types of flexion are comparable and are up 
to two orders of magnitude tighter than those from shear. Flexion therefore appears 
to be a very promising potential tool for constraining the shapes of galaxy haloes from 
future surveys. 
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1 INTRODUCTION 

Light from distant galaxies is perturbed by the curvature 
in space-time induced by intervening matter. This distorts 
the galaxy images and is referred to as gravitational lens- 
ing. In the most extreme cases of strong lensing the galaxy 
images are bent into thin arcs around a concentration of 
mass, and there may be multiple images of the object. 
Much more often a galaxy image is distorted only very 
slightly, and the image is sheared by a simple matrix dis- 
tortion. Thus an intrinsically circular object will appear 
as an ellipse. This first-order effect is referred to as weak 
gravitational lensing. The next order gravitational lens- 
ing effect, known as "flexion ", is a relatively new way of 
probing gravitational lenses ( Goldberg & Nataraian 2002; 
Irwin fc Sh makova 2003; Bac on et al. I l200rj ; Schneide r fc Erl 
20081 ; Bernstein fc Nakaiimal 20081 ). Other second-order ef- 
fects are not expected to be produced by simple gravita- 
tional lensing_j2ut_could_be a signature of systematic ef- 
fects ((Bacon fc Schaeferllioosl ). 

Flexion comprises of the one-flexion, or "displacement", 
and three-flexion, or "cardioid shift" . When combined with 
the shear they describe how "banana-like" a galaxy behind 
a gravitational lens appears. A disadvantage of the con- 
ventional weak lensing shear technique is that galaxies fre- 
quently have an intrinsically elliptical or stretched shape, 
and the extra stretching due to lensing is small by compar- 
ison. This effectively introduces extra noise into the mea- 
surements and requires an average over many galaxies to 
observe a significant lensing effect. However, galaxies are 
much less often intrinsically banana-shaped, and thus flexion 
measurements are expected to have a much higher signal to 
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noise ratio in some regimes. In addition the information ob- 
tained is orthogonal to that obtained from shear. Therefore 
information about lens systems derived from flexion can be 
combined with that from shear to obtain tighter constraints. 

Techniques have been developed to measure flexion 
from astronomical images, in the p r esence of pixelisation 



and noise dOkura et al.l 120071. 120081 
120071 ; iMassev et al.l 1200 
ion has been used to 
mass in clusters of galaxies 



Goldber g fc Leonard! 
Irwin fc Shmakovall2006|) . Flex- 
constrain the distribution of 



(llrwin fc Shmakoval 2003; 
iLeonard et ai1l2007l ; lOkura et al1l2008l ; ILeonard et al.ll2008l ) 

and gi ve indications of the distribution of mass in blank 
fields l|lrwin fc Shmakoval [200^ ; llrwin et~alll2007l ). Flexion 
due to the mass distribution in foreground galaxies (galaxy- 
galaxy flexion) has been observed in the Deep Lens Survey 
l|Goldberg fc Baconll2005l ). In this paper we discuss the use 
of galaxy-galaxy flexi on to measu r e halo ellipticities. This 
builds on the work of iBacon et alj l|2006l ). 

A method for using weak gravitational lensing shear 
to determine galaxy halo ellipticit ies by stacking many 
galaxi es together was proposed in iNataraian fc Refregierl 
(2000). This has been used to determine the ellipticities 
of ga la xy haloes llHoekstra et all |2004| ; iMandelbaum et al.l 



120061 : IParker et all 120071) 



haloes l|Evans fc Bridle! |200S| ) 



and cluster dark matter 
Ellipticities of individ- 



ual clusters have a lso been measured dCvpriano et al 1 I2004J : 
ICorless et alll2008h . 

Accurate determination of the shapes and orientations 
of galaxy dark matter haloes can provide constraints on 
models of galaxy formation and the nature of dark mat- 
ter. N-body simulations of non-interacting cold dark matter 
predict th at the these haloes s hould be triaxial prolate ellip- 
soids (e.g. lAllgood et al1l2006l ). Simulations that use cooling 
gas and dark matter p redict that dark matter haloes should 
be more spherical (e.g. iKazantzidis et al.ll2004l ). Predictions 
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in the Modified Newtonian Dynamics (MOND) paradigm 
predi ct that dark matter haloes should look isotropic at large 
radii (jMortlock fc Turnerl200lh . Previous weak lensing stud- 
ies have found flattened, elliptical galax y dark matter haloes 
which disfavor the MOND paradigm jHoekstra et alj|2004 
iMandelbaum et al.ll2006l ; iParker et al.ll2007T ). Thus through 
measuring the shapes of dark matter haloes, galaxy-galaxy 
lensing can provide constraints on galaxy formation models 
and the nature of dark matter. 

In Section[2]we derive expressions for the components of 
the flexion produced by an elliptical gravitational lens of ar- 
bitrary dimensionless surface density profile and illustrate 
the flex ion induced by an e lliptic al Navarro- Frenk- White 
(NFW; iNavarro et ail 1 19961 . 1 19971 ) ) lens. In Section |3] we 
show predictions for the uncertainties expected on the mean 
ellipticity of galaxy dark matter haloes. We conclude in Sec- 
tion n 

We adopt a fiducial flat cosmological constant domi- 
nated cold dark matter (ACDM) cosmology with Om = 0.3, 
SIa = 0.7 and a fluctuation amplitude parameterized by the 
root mean square density perturbation in 8 h~ x Mpc spheres 
at the present day in the linear regime, ag = 0.74. 



2 FLEXION BY ELLIPTICAL HALOES 

N-body simulations predict dark matter haloes to be tri- 
axial ellipsoids as opposed to spherical. A homeoidal ellip- 
soid is one in which shells are concentric triaxial ellipsoids 
all w ith the same ellipticity and orientation (e.g. ISchramml 
1990). Thus homoeoidal ellipsoids maintain their symme- 
try when projected onto a plane. In practice dark matter 
haloes are more complicated and may contain substructures 
and twisted isodensity contours and hence would not be cor- 
rectly described by a homoeoidal density profile. However, 
in galaxy-galaxy lensing many lens haloes are compiled on 
top of one another, ideally with their major axes aligned. 
This would smooth out any twisting or substructure, thus a 
homoeoidal description of the collated density profile is jus- 
tified. For a detailed disc ussion on various ellips oidal models 
and related theorems see IChan drasckhar ( 19691) . 

We now consider flexion by an elliptical mass distri- 
bution. An elliptical potential corresponds to an approxi- 
mately elliptical mass distribution in the small ellipticity 
limit. For large ellipticities t he m ass distribution becomes 
dumbell shaped. iBaltz et al.l (|2009l ) get around this by find- 
ing a particular sum of elliptical potentials that roughly 
mimics an elliptical mass distribution. Here we derive the 
exact equations for an elliptical mass distribution for the 
first time. 



2.1 First and Second Flexions 

The weak lensing shear and flexion are conveniently de- 
scribed ushig_a£cmplex_space formalism introduced to flex- 
ion by iBacon et al.l (|2006l ) that allows the magnitude and 
direction of the distortion effects to be expressed simulta- 
neously. The shear is composed of real and imaginary parts 
7 = 71 + 172 where 71 describes the stretching along the 
x and y axes, and 72 describes stretching along the diago- 
nals of the coordinate system. Similarly we can write for the 



one-flexion 

T = T\ + 1T2 — 9*7 
and for the three-flexion 
Q = Q\ + iQi = d-y. 



(1) 



(2) 



where d denotes a complex derivative o perator d = d/dx + 
id/dy with <9* as its complex conjugate (IBacon et al.ll200"6l ). 

The shear is a second derivative of the lensing potential 
ip; 71 = ^(tpxx+tpyy) and 72 = ipxy, where subscripts denote 
partial differentiation. The real and imaginary components 
of the flexion can therefore be expressed explicitly in terms 
of the third derivatives of the potential 



(3) 
(4) 
(5) 
(6) 



So in order to calculate the flexion components induced at 
position (x, y) on the image plane by an elliptical halo all 
one has to do is find derivatives of the potential around an 

elliptical lens. 

By differentiating equations 9-11 of iKeetonl l|200ll ) we 
obtain expressions for the third derivatives of the deflection 
potential (equations [7] - I10[) induced by an elliptical dark 
matter halo of arbitrary lensing convergence, k. In general 
relativity and the thin-lens approximation, the lensing con- 
vergence is simply the mass density projected along the line 
of sight. Here we are mapping a one-dimensional k profile, 
onto homoeoids with minor to major axis ratio q. We find 



Tx = 




T 2 = 


— (ipxxy + Ipyyyji 


Qi = 


— {ipxXX — "Slpxyy j 


Q2 = 


— (Slpxxy — 'tpyyyj 



Tpxxx — 

Yvvy — 

Ipxxy — 

tpyyx = 

Here 

K n (x,y) 
and 



6qxKo + 4qx Lo, 
6qyK 2 + 4<?2/ 3 L 3 , 
2qyK\ + 4qx 2 yLi, 
2qxK\ + 4qy 2 xL 2 



1 uk! (£(«)) 2 du 

[1 - (1 - Q 2 ) u ]n+l/2 



Ln(x,y) = 



1 u 2 k" {j{u)fdu 
[1 - (1 - g 2 )u]«+V2 



(7) 
(8) 
(9) 
(10) 



(11) 



(12) 



are one-dimensional integrals from the centre to the ellipse 
intersecting the (x, y) position of interest. These are given in 
terms of the first and second derivatives of the convergence 
as a function of ellipse coordinate 



where 



[1- (i-g 2 H 



(13) 



(14) 



is the elliptical coordinate. 

By substituting equations [71 - HOl into [3l to l6l we can ob- 
tain the first and second flexions at any point on the plane 
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Figure 1. Illustrations of the three simplest lensing distortions for an NFW dark matter halo with an ellipticity of e = 0.3 and mass of 
10 12 h^ 1 Mq and concentration parameter c = 7.2. The top left panel shows arrows pointing in the direction of the one-flexion vector, 
the length of the arrow is proportional to the magnitude of the one-flexion. The top middle panel is a colour map showing how Tr varies 
across the image the top right panel shows how Tb varies. The three panels in the middle row illustrate the equivalent quantities for 
the shear and the bottom row for three-flexion. Colour bars on the right of each flexion image give the distortion per arcsecond. The 
contours on the left hand panels are at convergence, k values of 0.01 and 0.02. Arrows and sticks that would be distractingly large have 
been removed. 



(x, y) arou nd an e ll iptica l lens for any k. Both ISchramml 
(1990) and iKeetonl 0001) note that only very simple mat- 
ter distributions give simple analytical relations for ellipti- 
cal lenses. However, f or realistic mat t er dis t ributions, such 
as the NFW profile (|Navarro et all (|l996h iNavarro et all 
{ 19971 )1; the expected flexions around a given lens can easily 
be computed. 

2.2 Radial Flexion 

The shear around an elliptical lens can be described in terms 
of the tangential 7t and cross 7x components 

7T = | 7 |cos(20;), (15) 
lx = | 7 |sin(20;) (16) 

where 

6?; =6? 7 (17) 

and 71 = |7|cos(2# 7 ), 72 = I7I sin(2# 7 ) and x = rcos(<^), 
y = r s'm(4>) where r 2 = x 2 + y 2 . Therefore 9 1 is the angle of 



the shear stretch to the coordinate axes x, y, f? 7 is the angle 
of the shear stretch relative to a rotated set of axes aligned 
along the radial direction and <fi is the angle subtended at 
the center of symmetry from the positive x axis to the vector 
x,y to the position of the background galaxy. 

We propose that the equivalent quantities for the flex- 
ions should be defined 



Tr = |^|cos(6»^), (18) 

Tb = |^"|sin(^), (19) 

Qr = |£|cos(30g), (20) 

Qb = |e|sin(30g), (21) 

where 

9'r = - + 7T, (22) 

O'g = Og-fr (23) 



and Tx = \F\ cos{$r), = \T\sm(Q T ), Qx = \Q\ cos(36» g ) 
and Q2 = 1 1/| sin (36*g ) . Here O'-p and 8'g are defined in such 
a way as to make Tb, and Qr positive around a circular 
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Figure 2. Number of galaxies per unit redshift as a function of 
redshift for a fiducial space-based survey. The green line indicates 
the cut-off between the lens and source population, shaded red 
and blue respectively. The blue dotted line shows the median 
redshift of the lens population whilst the magenta dotted line 
shows the median redshift of the source population. 



lens. These 'radial' components are the gravitational contri- 
butions of the flexion around a circular lens, any Tb or Qb 
around a circular lens is hence bogus. 

We illustrate Tr, Tb, 7t, 7x, Qr an d Qb for an el- 
liptical NFW profile in Figure [1] For an elliptical singu- 
lar i sothermal sphere prof ile the shear is entirely tangen- 
tial l|Kormann et al.lll994 ). just as for a circular lens. We 
see from the central panel that the result for an NFW pro- 
file has a small but non-zero cross shear. The cross shear is 
such that the direction of the shear sticks at positions on the 
y — x line are slightly more horizontal than a purely tangen- 
tial line. They therefore tend towards alignment along the 
isodensity contours. 

The one-flexion has a significant cross component. This 
is clearly seen in the top left stick plot, in which the arrows 
have a vertical component and therefore are more closely 
perpendicular to the isodensity contours. The three flexion 
has a small cross component and the far leg of the trian- 
gle is close to radial despite the ellitpicity of the lens. The 
exact perturbation from radial varies according to the angu- 
lar position on the lens, but is mostly pointing even further 
from the radial line than the perpendicular to the isodensity 
contours. 



3 APPLICATION TO SURVEYS 

Here we estimate the potential for future surveys to mea- 
sure the ellipticity of haloes using stacked shear and flexion 
maps. We calculate flexions and shears on a grid covering 
one arcminute square centred on the lens centre. This is 
analagous to stacking galaxies in postage stamps of one ar- 
cminute square, centred on the halo center with the major 
axes aligned along the x axis and averaging observed flex- 
ions, or shears, in pixels. This method is most applicable to 
large scale weak lensing surveys. 



3.1 Lens and source populations 

We need an estimate of the lens parameters, and of lens and 
source galaxy redshifts, in order to calculate the flexions 
we would expect to observe at positions around the lens. In 
this work we carry out an illustrative calculation using mean 
quantities for each of the lens and source populations. 

To estimate the typical redshift of the lens and source 
we divid e up the redshift distri bution of galaxies from Equa- 
tion 6 of iBlake fc Bridle! 1)20051 1 



dN 
dz 



^ z i 
oho — o exp 



(I)' 



where the characteristic redshift 



z = 0.055(r ; . lm - 24) + 0.39 
and overall surface density 
_ 35400 ,r Um , 



(24) 



(25) 



(26) 



are found from linear fits to their Table 1. Note that this 
table is derived from observed luminosity functions and ex- 
tends only up to a magnitude of 24, therefore the above 
equations are used to extrapolate beyond the numbers pro- 
vided. For general ground-based surveys we use a = f3 — |. 

We define the lens population such that there is one lens 
per square arcminute postage stamp. The postage stamp size 
would ideally be small to avoid overlaps between lenses, but 
to get a good signal from a finite survey area we make a 
compromise at 1 arcminute. Beyond this the flexion signal 
is very noisy. We assume that the lens population of inter- 
est will be all the low-redshift galaxies, up to some redshift 
cut-off. We take the cut redshift to be such that there is 
one lens galaxy per square arcminute. This means that in 
our estimation the number of lenses in any survey is solely 
dependent on the area of that survey whereas the number 
of sources is dependent upon both the depth and area of the 
survey. The redshift of the background galaxies is taken to 
be the median redshift of the population posterior to this 
cut-off. Similarly the lens redshift is taken to be the median 
redshift of the lens population. 

Using the median redshift of the lens population we 
can estimate the lower mass limit of the observable haloes 
at this redshift using the Sheth-Tormen mass function 
jSheth fc Tormenlll999 



n h (M)dM = A 



l + — q 



2 p m dv 
7t M ~d~M 



exp 



y W(27) 



where v = ^/aSrD(z)a{M), A = 0.322, a = 0.707 and 
q — 0.3. We do this by varying the lower mass limit until 
the number density of haloes matches the number density 
of galaxies for our given magnitude limit. This mass is then 
used as our lens mas, M200, in our calculations of the lens- 
ing potential. Technically this is the lower limit to the mass 
rather than the mean mass. However it will be close to the 
mean due to the shape of the mass function. The small dif- 
ference makes our uncertainty estimates on the conservative 
side. 

Hypothetically any k profile could be plugged into[7]-[TD] 
but we have chosen an NFW profile because it is reasonably 
realistic. The value of the NFW concentration parameter is 
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0.32 



Ellipticity Ellipticity Ellipticity 

Figure 3. Probability as a function of lens ellipticity for the fiducial large area space-based survey. The true ellipticity is shown by the 
red dashed line and the blue dotted lines show the 68% confidence limits. The left panel shows the constraints from one-flexion, the 
central panel for three-flexion and the right hand panel for shear. 



estimated by assuming that its relati onship with A /200 has 
the functional form of equation 12 of ISeliakl (|2000h 

( M 20 o \ B 

c = c \-m7) ' (28) 

where Co = 10, B = —0.15 and M* is a function of z. 

As magnitude limit is increased, more low luminosity, 
low mass, low redshift galaxies become visible. If these are 
included in the lens population then the median redshift 
and thus assumed mass of the lens galaxy population will 
decrease. This decreases the lensing efficiency resulting in 
little improvement on potential dark matter halo ellipticity 
constraints as the depth of the survey is increased. Therefore 
by default we carry out calculations in which the lens pop- 
ulation only contains galaxies brighter than a magnitude 
limit of 24, even though the survey itself may be deeper. 
This means that the properties of the lens population are 
unchanged as the survey depth is increased. This includes 
the lens cut-off redshift. The source population is then taken 
to be all the galaxies at redshifts above the lens cut-off red- 
shift and down to the magnitude limit of the survey. The 
source population redshift then increases with magnitude 
limit giving a stronger flexion signal. 

The lens population from an ru m — 24 survey covering 
20,000 square degrees of sky has a median redshift of zi = 
0.19 and there are n s = 6 x 10 8 sources with a median 
redshift of z s = 0.58. We find A/200 = 2.6 x IO^A/q/i -1 
with concentration parameter c = 22. An ru m = 26 survey 
has n s — 3.1 x 10 9 sources with z s — 0.72. Whereas if we 
were to allow the lens population to change as the survey got 
deeper, we would have in the ru m = 26 survey a lens redshift 
of zi = 0.14, a lens mass of A/200 = 9.3 x 1O 9 A/ /i" 1 with 
a concentration parameter of c = 26. There would be n s = 
3.2 x 10 9 sources with a median redshift of z s — 0.71. These 
subtle changes add up to a noticeable change in ellipticity 
constraints (Figure 

We also show results for a fiducial large-area space- 
based survey covering 20,000 square degrees of sky with 35 
galaxies per square arcminute. We assume a = 2, f3 = |, 
20 = 0.63 and So = 27. This has lenses with a median red- 
shift of zi — 0.19 and n s = 2.5 x 10 9 sources with a median 
redshift of z s = 0.91. This gives M 20 o = 2.6 x IO^Mq/i" 1 
with c = 22. 



3.2 Constraints on Lens Ellipticity 

We now use the lens and source redshifts and number densi- 
ties described in the previous subsection to make simulated 
stacked shear and flexion maps with uncertainties. We fit 
unknown lens parameters to these maps to obtain probabil- 
ity distributions as a function of lens ellipticity. We compare 
the results from shear, one-flexion and three-flexion. We also 
investigate how our results depend on the survey depth. 

The final input to our analysis is an estimate of the un- 
certainty on the shear and flexion measurements for a single 
background galaxy. We assume the uncertainty on the shear 
measurement for each galaxy is 0.3 for each shear compo- 
nent. Although this number does depend on the depth of 
the survey, for deep surveys it is dominated by the intrin- 
sic distribution of galaxy ellipticities which contribute about 
0.2 to the shears. Owing to the extreme di ficulty in obt ain- 
ing flexion measurements from real data ([Rowc 20q3) we 
have used a conservative estimate of ot = og = 0.1. This is 
larger than the estimat ed noise on the flexions <jt = 0.03 
and og = 0.04 used in lGoldberg fc Leonard] (120071 ). Again 
it is known that noise on the flexion is dependent on survey 
depth. These uncertainties are divided by the square root of 
the number of galaxies in each pixel on the collated postage 
stamp to obtain an estimate of the noise for that pixel. 

We consider trial values of the lens mass and lens el- 
lipticity and compute predicted shear and flexion maps. We 
keep the lens centered on the origin and the orientation along 
the i-axis since this would be implicit in the stacking proce- 
dure. We compute the probability of each model by calculat- 
ing a x 2 between the trial map and the fiducial map, using 
the uncertainties computed for each pixel. We marginalise 
over the mass to find the probability as a function of lens el- 
lipticity and find the uncertainty on the ellipticity from the 
68 per cent confidence limits. This is illustrated in Fig. [3] 
for the fiducial space-based survey. The constraints on el- 
lipticity from the space based survey are o e = 2.0 x 10 -5 
from one-flexion, <r e = 9.1 x 10 -6 from three-flexion and 
o e = 4.5 x 10 -3 from shear. It is clear that flexion is a much 
more powerful method with which to constrain lens ellip- 
ticity than shear alone. Indeed it is not even convenient to 
show the constraints from shear on the same axes. 

We show the uncertainties in the lens ellipticity using 
flexion estimated for ground-based surveys as a function of 
survey depth in Fig. [4] We see that the two types of flex- 
ion give similar size constraints, with the three flexion giving 
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0.5 1 1 1 1 

24 24.5 25 25.5 26 

Figure 4. Error on cllipticity with respect to magnitude limit of a 
ground based survey. The dotted lines show constraints obtained 
using a varying lens population whilst the blue solid line and 
red dashed line show constraints using the lens population from 
the ru m = 24 ground based survey. As the magnitude limit of 
the survey increases so the uncertainty in the mean ellipticity of 
galaxy dark matter haloes decreases. There is about a factor of 
two improvement in moving from an r; im = 24 to an ri im = 26 
survey. 




ajr/ag/arcsec 

Figure 5. Error on dark matter halo ellipticity as a function of 
error on flexion measurements for one-flexion (dotted line) and 
three-flexion (solid line). This is for our space based survey. The 
estimated error on the flexions o f = 0.03 and ag = 0.04 
con sistent with iBacon et al.l ||2006| 1 . [Goldberg fc Leonard! (120071 1 
and iRowd (|2007F are illustrated, as is the value we have used of 
Ujr = ag = 0.1. 



slightly tighter constraints. When the lens population is kept 
fixed as the survey depth is increased (solid lines) the uncer- 
tainty approximately halves as the magnitude is increased by 
two. When the lens population is varied to keep the number 
density of lenses constant the constraints vary little. This is 
because, as discussed earlier, the lenses become closer and 
lower mass. ,j,From our rj jm = 24 survey a e — 4.7 x 1CP 5 
with one-flexion and a e — 4.4 x 10~ 5 with three flexion. 
This compares to a e — 8.2 x 1CP 3 using the shear. From 
our rum = 26 survey constraints are a e = 2.0 x 10~ J with 
one flexion and a e — 9.7 x 10 -6 with three flexion. Shear 
provides constraints of a e — 4.4 x 10 -3 . These numbers are 
very similar to those from the space-based survey. 

In Fig. 13.21 we illustrate the ellipticity uncertainty as 
a function of assumed flexion uncertainty. This shows that 
the relative error on the flexion measurements themselves 
can be large without compromising flexion's ability to con- 
strain the ellipticity of dark matter haloes. It has been ob- 
served that at fainter magnitudes the uncertainty in flexion 
measurements is greater, this effect is most noticeable for 
three-flexion (David Bacon private communication). If, for 
example, the uncertainty in one-flexion measurements were 
~ 0.1 and in three-flexion ~ 0.3 then the constraints ob- 
tained from each would be roughly the same. 



4 CONCLUSION 

We have generalised the equations for shear from an el- 
lipsoidal mass distribution to flexion. These could in prin- 
ciple be extended to even higher order lensing effects 
by further differentiation. This differs from previous work 



bv lBaltz et al.l l|2009l ) who approximate an ellipsoidal mass 
distribution with a sum of ellipsoidal potentials. 

We have demonstrated that galaxy-galaxy flexion can 
constrain the ellipticity of galaxy dark matter haloes up to 
two orders of magnitude better than shear. Such precision 
measurements would be very useful in constraining models 
of galaxy formation and galaxy dark matter. Increasing the 
magnitude limit of a ground based survey from ru m = 24 to 
rum = 26 tightens the constraints on ellipticity by up to a 
factor of two. 

The exact values depend heavily on the assumed in- 
trinsic flexion uncertainty coming from the intrinsic galaxy 
shapes. As mentioned in section [3~2l the flexion uncertainty 
is expected to increase with survey depth, this could be ac- 
counted for with a more detailed treatment. This will make 
our predicitions for ellipticity constraints slightly optimistic. 
Furthermore, we have not considered the CCD resolution for 
the space based survey. Pixels much larger than an HST style 
camera could substantially increase flexion uncertainty. Ad- 
ditionally, we have assumed Gaussian errors on the flexion 
measurements yet it is known tha t observed flexion errors 
are far from Gaussian (|Rowell2007| ). 

In carrying out a galaxy halo ellipticity observation 
there is considerable freedom in selecting the lens and 
source population to include in the analysis. Since we are 
interested in a raw comparison between shear and flex- 
ion we used a simple redshift cut on the lens population. 
However it may be useful to split t he lens population by 
observed ellipticity before stacking (jHoekstra et al.l |2004| ; 
iMandelbaum et alJlioOrj ; IParker et al. 1 120071 1. to increase the 
signal-to-noise of any detection of a difference between the 
halo and light ellipticities. Furthermore it is of interest 
to study the evolution of halo ellipticity with halo mass 
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and redshift 



Hopk ins et al.l 



Kasun fc Evrardl 120051 ; lAllgood et all 120061: 
20051 : iHo et alj|2006h . 



In practice the selection of lens and/or source galaxies 
may be carried out using photometric redshifts, which would 
not allow a straight cut in redshift, but would cause the two 
populations to overlap somewhat. This must be carefully ac- 
counted for in determinations of the mass, but has a second 
order effect on the ellipticity. 

We have further assumed that galaxies at all positions 
behind the lenses may be used for shear and flexion anal- 
ysis. However it will be impossible to measure background 
galaxy shapes very close to the foreground galaxy, due con- 
tamination of the measurement by light from the foreground 
object. Depending on the ellipticity and orientation of the 
light relative to the mass, this will impose a roughly ellip- 
tical mask on the background galaxy catalogue which will 
weaken both the constraints on the halo ellipticity for both 
shear and flexion measurements. The exact extent will de- 
pend on the resolution of the telescope and the efficiency of 
the shear or flexion measurement code to deal with varying 
background levels across the image of a lensed galaxy, which 
is not yet widely established. 

Current surveys such as the Hubble Space Tele- 
scope (HST) Cosmic Evolution Survey (COSMOS) and the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) 
might already be able to place tight constraints on halo el- 
lipticity using flexion, if the flexions can be measured with 
sufficient accuracy. We look forward to future measurements 
from high precision imaging surveys such as the Dark En- 
ergy Survey (DES), Pan-STARRS, the Large Synoptic Sur- 
vey Telescope (LSST) and a space mission such as Euclid 
and/or the Joint Dark Energy Mission (JDEM). 
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